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The three main contributions to the nuclear Hamiltonian — monopole, quadrupole and pairing — 
are analyzed in a shell model context. The first has to be treated phenomenologically, while the other 
two can be reliably extracted from the realistic interactions. Due to simple scaling properties, the 
realistic quadrupole and pairing interactions eliminate the tendency to collapse of their conventional 
counterparts, while retaining their basic simplicity. 
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Nuclei are systems of interacting particles. Everybody 
agrees with this statement. But then: which interaction, 
which particles? These questions raise some difficulties, 
so the first statement is replaced by : Nuclei are systems 
of quasi-particles, interacting via effective Hamiltonians 
(or Lagrangians) . The questions recur, and now the an- 
swers vary widely. Restricting attention to theorists who 
are content with the idea that the particles are basically 
neutrons and protons, there main possibilities concerning 
the interaction: 

R. The realistic approach consists in extracting it from 
NN data, and then take seriously the idea that the many- 
body Schodinger has to be solved exactly. 

P. The phenomenological approach (es) derive the in- 
teraction from the nuclear data they are supposed to ex- 
plain. There are two variants: 

MFP Mean field phenomenology restricts to a maxi- 
mumn the number of parameters, and directs atten- 
tion to global properties and general spectroscopic 
trends. 

SMP Shell model phenomenology is prepared to intro- 
duce as many parameters as necessary to explain 
spectroscopic detail. 

Conceptually, R is the most satisfactory but, in ad- 
dition to its difficulty, it has suffered so far from a ma- 
jor drawback: The realistic interactions do not saturate 
well, i. e., they do not give the right binding at the 
right radius. Furthermore, they fail to produce the cor- 
rect shell structure: doubly magic closures are missed. 
Much progress has been made recently in fitting the NN 
data perfectly, and in calculating exactly for very light 
nuclei. However, the problems remain, and phenomeno- 
logical three-body forces have to be called in to solve 
them. 

Historically, P, has enormous importance, and both its 
variants rest on pioneering work that remains of perma- 
nent value. It can be summarized by saying that, what- 
ever their guise, the effective Hamiltonians must contain 
a monopole term that produces the spherical closures, 
a quadrupole one that induces deformation, and a pair- 
ing force that favours condensation at the vicinity of the 
Fermi surface. 



My purpose is to show how these things are possible, 
within a framework that starts from R. The idea is that 
the trouble in realistic interactions is concentrated in the 
monopole part, which can be separated rigorously from 
the rest. Moreover, it is formally very simple and can be 
extracted from the data with a mimimum of parameters. 
Once this is done, the rest of the Hamiltonian (which 
we call multipole), turns out to be independent of the 
type of realistic interaction used. What's more, it does a 
magnificent job in large scale shell model calculations. 

The aim of these notes m is to provide a unified in- 
troduction to references [y,^ , dealing with the monopole 
and multipole Hamiltonians respectively. 

NOTATIONS A few equations have to exhibit explic- 
itly angular momentum (J), and isospin (T) conserva- 
tion. I will use Bruce French's product notation B) 

F stands for JT. Then {-f = (-)"'+^, [r] = 
(2J-I- 1)(2T-|- 1), and in general F(r) = F{J)F{T). Or- 
bits are called r, s, etc., and I use (— )'' = (— )J'-+i/2^ 
[r] = 2{2jr + !)• Expressions carry to neutron-proton 
formalism simply by dropping the issospin factor 

rrir is the number of particles in orbit r, T^ is used for 
both the isospin and the isospin operator. In neutron- 
proton {np) scheme, rurx specifies the fluid x. 

^rsV is ^^ operator of type a\a\ coupled to good spin 
and isospin JT. S^^ is an operator of type ajog coupled 
to good spin and isospin At. 

l/^i„ is a two body matrix element. W^^^^ is used after 
the monopole part has been subtracted. 

p is the principal oscillator quantum number. 



I. THE MONOPOLE HAMILTONIAN 

The only assumption that will be made is the existence 
of an effective potential smooth enough to do Hartee Fock 
(HF) variation, and capable of yielding good results. 

Now: given a Hamiltonian 7i, (/C is the kinetic energy) 
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it is always possible to extract from it a monopole part 
TC^, whose expectation value for any state is the average 



energy of the configuration to which it belongs (a config- 
uration is a set of states with fixed rrirx for each orbit). 
In particular, WJ^reproduces the exact energy of closed 
shells {cs) and single particle (or hole) states built on 
them ((cs) ± 1), since for this set [cs ± 1) each configura- 
tion contains a single member. The result is standard 0, 
and we simply write it, {d stands for diagonal) 



'Ht =^"+ Y. Vr^'mU^sx' - SrsSxx' 



(2) 



which reproduces the average energies of configurations 
at fixed mrx'mrx' ■ and isospin (jnT) 
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which reproduces the average energies of configurations 
at fixed rrirTr. 

Using Dr — 2jr + 1, we rewrite the relevant centroids 
incorporating explicitly the Pauli restrictions 
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In the np scheme each orbit r goes into two rx and rx' 
and the centroids can be obtained through {x ^ x') 
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Under Hartree Fock variation, the rrir and Tr operators 
will go into non diagonal ones of the type S^^. Therefore 
TiJ^^and TiJ^j^should be generalized to the full monopole 
Hamiltonian, Tim, containing all two body quadratic 
forms in the S^[ operators, see [p|,0). The task is not 
trivial, and we can avoid doing the variation explicitly, 
as we see next. 



ViLu)khnn - —ViuJo)klrnn, hw = 3A.6A^/y{r^, (8) 
LUq 

where the second equality is adapted from |8[| . A very ac- 
curate fit to the spherical radii (a variant of those in g]) 
yields (r^) = 0.89^^^/3 ^j^h 

p = Ai/3(l-(2T/A)2)e(3-5M), so hw^iQ/p. (9) 

This is very simple, but it cannot be the full 
story, because the competition between potential 
and kinetic energy (which has the same scaling, 
/C* = ilT = fixj/2^ (p-t- 3/2)mp, irip is the number of 
particles in harmonic oscillator (HO) shell p) would lead 
to a trivial equilibrium at fiij = cx) or 0. The terms that 
cannot scale strictly with lu must be those that go as 
the total number of particles A, since they are alone re- 
sponsible for saturation. We shall not try to discover the 
saturation mechanism, but simply note that Hf^should 
produce (the right) smooth contributions in A and A^/"^, 
and a symmetry energy in T{T+l)/A and T{T+l)/A'^^^. 
Next we separate these "liquid drop" terms, from those 
that produce shell effects Hm—Tif'ri^ + Ti^, and set out 
to parametrize efficiently Ti^. 



1. Tifn o,nd the collective monoplole Hamiltonian W 

By definition, shell effects go as A^^"^, and we have as- 
sumed they can be separated those in A and A^^^ . Struti- 
nsky's famous theorem tells us this is possible. Here I 
will show an explicit example, which is the starting point 
in the construction of 7i^. Independently of the detailed 
saturation mechanism, there must be "something" in 0.^ 
capable of cancelling K. As it is always possible to sepa- 
rate from Hm a leading term in m{m—l)/2 = A{A—l)/2, 
the "something" must be related to it. As the end result 
must go like A, this term must be scaled by A~^. How- 
ever, there is no such thing in 0.^ if wc take the scaling to 
he huj , i. e., A~^/^. By using the techniques of section ||, 
we discover that the term we are after is 



M^ = ^(mp/, 



Dr, 



(p+l)(p + 2), (10) 



A. Scaling 

HF variation is necessary to ensure that the system as- 
sume its correct radius, at an energy close to the correct 
one. In all representations 7if„has the same form, that 
must change smoothly as N and Z change. To discover 
how the evolution takes place, we first note that matrix 
elements for a potential of short — but non zero — range 
scale as the oscillator frequency huj (a consequence of the 
Talmi Moshinsky transformation). Then we can write 



which has the remarkable property that huj{W — 4:K)/4 
cancels exactly to order A^/"^ and (as we shall see) pro- 
duces strong shell effects. W is called the collective 
monopole Hamiltonian for reasons that will become fully 
apparent in section II A. 

W — AK is taken to be the first contibution to 7if„, 
which has been fitted on the set of known (cs)±l spectra. 
For the explicit construction of the rest I refer to . For 
the present purpose it will be sufficient to mention that it 
includes axil- s + 1-1 piece that behaves as one-body and 



has the same coefBcient as W — AK. Then, there are two- 
body terms that ensure the correct evolution of {cs) ± 1 
spacings from HO to EI closures (extruder-intruder, such 
as N Z = 14, 28, 50, etc). Fig |i| illustrates the sheU for- 
mation mechanism for T — b nuclei. Note the beautiful 
HO closures generated hyW — iK, erased hy I- s + l-l (as 
in ^°Ne at iV = 20), and replaced by EI closures through 
the two body terms. 



Replacing pairs by single indeces rs = x, tu = y in 
eq. (n2h and rt ^ a, su = b in eq. (O), we bring the 
matrices W^y and /2(, = frtsu = CrsCtu, to diagonal form 
through unitary transformations U^ 



U-^WU = E ■ 



W. 



xy 



k' "-ak- 



z2^^kUykEl 



(16) 




15 20 25 30 35 40 45 50 55 60 
FIG. 1. Shell formation mechanism in T = 5 nuclei. Nega- 
tive binding energies in y-axis are repulsive. N in x-axis 



II. THE MULTIPOLE HAMILTONIAN 

The multipole Hamiltonian is defined as Ti.M—'H-Ti.m- 
As we are no longer interested in the full H, but its 
restriction to a finite space, T^Mwill be more modestly 
called Hm, with monopole-free matrix elements given by 



W, 



JT 



K 



JT 



^rsStiiV^. 



(11) 



To concentrate on the physics I shall only write the for- 
mulae that are strictly necessary to show how the pairing 
plus quadrupolc interactions appear in their normalized 
form of. Abundant detail is given in ||3|. 

There are two standard ways of writing Hm'- 



Hm — 



E 



wt 



i^rsT • ^tuF, or 
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Hm^ 5] [7]^/V,,J57t5L)°, 
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where Q^s = (1 + ^rsf'^'^ 1^ j and the matrix elements are 
related through 

^?...^EH— {::n<.[r], (14) 



w^.^.*.-EH^^*"^-^{;;; !j}--« 



7]. (15) 



and then, 

fc,r X y 



H 



M 



E' 

/c,7 



,7 C7 



b 



^bk'-'b 



h] 



1/2 



(17) 



(18) 
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which we call the E and e representations. Since Tim 
contains all the 7 = 00 and 01 terms, for Hm, uJr^u — 
oJrstu = 0- There no one body contractions in the e rep- 
resentation because they are all proportional to oJrstu- 

The eigensolutions in eqs. (|l^) and (^ using the KLS 
interaction ||l^,|ll| — in ref. it is explained in detail why 
this venerable choice is a good one — for spaces of one and 
two major oscillator shells. The density of eigenvalues 
(their number in a given interval) in the E representa- 
tion is shown in fig. H for a typical two-shell case. 
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FIG. 2. _E-eigenvaIue density for the KLS interaction in the 
pf-|-sdg major shells hui = 9. Each eigenvalue has multiplicity 
[r]. The largest ones are shown by arrows. 



It is skewed, with a tail at negative energies which is 
what we expect from an attractive interaction. 
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FIG. 3. e-eigenvalue density for the KLS interaction in the 
pf+sdg major shells. Each eigenvalue has multiplicity [7]. 
The largest ones are shown by arrows. 
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FIG. 4. _E-eigenvalue density for the KLS interaction in the 
pf+sdg major shells huj = 9, after removal of the five largest 
multipole contributions. Each eigenvalue has multiplicity [F] . 
The largest ones are shown by arrows. 



The e eigenvalues have a number of simple proper- 
ties demonstrated in p, Appendix B] : their mean value 
always vanishes, their width is \/T/8 of that of the 
E distribution, and they are twice as numerous. In 



fig. 1^ we find that they are very symmetrically dis- 
tributed around a narrow central group, but few of 
them are neatly detached. The strongest have 7^^ — 
1-0, 1+1, 2+0, 3-0, 4+0. If the corresponding 
eigenvectors are eliminated from H in eq. ( pj| ) and 
the associated H in eq. (Iq) is recalculated, the E distri- 
bution becomes quite symmetric. (The boldface is meant 
to call attention to the one very bad erratum in H). The 
result is shown in Fig. H. The residual skewness is en- 
tirely accounted for by the F — 1+0, 0+1 and 2+0 peaks, 
whose strength is somewhat eroded but remains substan- 
tial. 

This result is most telling because from the work of 
Mon and French [|2| (see also the references quoted in Q| ) 
we know that a symmetric E distribution will lead to 
spectra in the m-particle systems that are identical to 
those of a random matrix. Therefore, we have found 
that - with the exception of three F peaks - the very 
few dominant terms in the e-distribution are responsible 
for deviations from random behaviour in Hm . Positively 
stated, these terms are at the origin of collective proper- 
ties. 

If the diagonalizations are restricted to one major shell, 
negative parity peaks are absent, but for the positive par- 
ity ones the results are practically identical to those of 
Figs. H and || except that the energies are halved. This 
point is crucial: 

// Upj and Up2 are the eigenvectors obtained in shells 
Pi and p2, their eigenvalues are approximately equal 

When diagonalizing inpi+p2, t/ie unnormalized eigen- 
vector turns out to be Up-^ + Up^ with eigenvalue e. 

In the figures the eigenvalues for the two shell case 
are doubled, because they are associated with normal- 
ized eigenvectors. 

To make a long story short: The contribution to Hm 
associated to the F = 01, and 7 = 20, 
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\{qp + qp+i) ■ (qp + qp+i), 



(20) 
(21) 



turn out to be (naturally) the usual pairing plus 
quadrupole Hamiltonians, except that the operators for 
each major shell of principal quantum number p are af- 
fected by a normalization. E^^ and e^" are the one shell 
values called generically e in the discussion above. To be 
precise 



-Tp — ^ ^ Zyj,^gj^SS^ /Sip , 


(22) 
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9p = E ^'s'Zra/A/'p, 
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where 




0, = jr + 1/2, qrs = V^dlr^r^lls), and 





f7p = ^a--2^: 






^Qr 



^(p+3/2)^ 



(24) 



More on the other collective terms of Hm in section III 



A. Collapse avoided 

The pairing plus quadrupole model has a long and glo- 
rious history [n3 14 1, and one big problem: as more shells 
are added to a space, the energy grows, eventually lead- 
ing to collapse. The only solution is to stay within lim- 
ited spaces, but then the coupling constants have to be 
readjusted on a case by case basis. The normalized ver- 
sions of the operators presented above are affected by 
universal coupling constants that dot change with the 
number of shells. Knowing that huio — 9 MeV, they are 
\E°'^\/huJa = g' = 0.32 and \e^°\/fiujo = k' = 0.216 in 
Eqs. pO^ and pl^. (One should not forget that for practi- 
cal use, these numbers must be renormalized. See H for 
details.) 

Introducing Amf ~ |(P/ + 3/2)'^, the total number of 
particles at the middle of the Fermi shell p/, the rela- 
tionship between g', k', and their conventional counter- 
parts M] is, for one shell 



0.32fta; 
0.216/10; 



19.51 

41/3/12/3 

1 216 
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= ^ = ^^-5/3 
2 ~ 2 



(25) 



What I propose next, is to show that both variants 
do indeed produce shell effects in A^/^, a necessity em- 
phasized at the beginning of section [Al. Then, by al- 
lowing the particles to be promoted to shells well above 
the Fermi one, it will become clear why the conventional 
forms produce collapse, and the new ones do not. 

Assume m = 0{Df) particles in shell pf, for which, 



D 



f 



0{A 



2/3N 



Consider first the pairing force, 
eracy D, it produces an energy 



On a space of degen- 



\G\ , 
Ep = m[D — m 



2) = -\G\0{mD). 



The first equality is a standard result. For D = Df, the 
conventional choice G — 0{A^^) leads to Ep = 0{A^/^), 
i. e., it guarantees the correct scaling for shell effects. 
That the result is not as trivial as it looks can be gath- 
ered from this quotation concerning the G coupling: "We 
know of no reliable way of predicting this A^^ depen- 
dence..." |3|. 

For a quadrupole force, an estimate for the en- 
ergy can be obtained by constructing a determi- 
nantal state that maximizes the quadrupole moment 



Qo = I]r=i(2^zi - ^xi - riyi), where n„,nyi,nzi are the 
number of quanta. The largest term in the sum is 
then 2p, the next 2p — 3, then 2p — 6, etc. Therefore 
Qo = 0{mp), and 



Er, 



-W\Ql^-W\0{m^D), 



which in turn explains (for D = Df) the origin of the 
usual choice x' = 0{A^^/'^) for the quadrupole strength, 
that leads to Eg = 0{A^/^). 

It is clear from equations (|2^) that the operators are 
affected by coefficients that go as A~^^^D~^ (instead of 
A'^) for pairing, and as A~^/^D~'^ (instead of ^~^/^) 
for quadrupole. For D — Df, the energies are again 
0{A^/^), but now this important empirical fact is a di- 
rect consequence of the interaction. For arbitrary D, the 
energies of the traditional (old) versions transform into 
the normalized (new) one as 



Ep{old) = O(^) 



^g(old) = 0( 



7n?D 



Ep (new) = 0{ 



Eginew) = 0( 



m , 



A^/^D' 



If the m particles are promoted to some higher shell 
with p = Pf + M, D « {pf + M)^, both energies grow 
in the old version. For sufficiently large M - because 
of the term in M^ - the gain will become larger than 
the monopole loss 0{mMhuj) ~ 0{MA^/^) that is only 
linear in M. Therefore the traditional forces lead the 
system to collapse. In the new form there is no collapse: 
Ep stays constant, Eg decreases and the monopole term 
provides the restoring force that guarantees that particles 
will remain predominantly in the Fermi shell. 

Let us go back now to the monopole term W that 
plays such a crucial role in TCf,^ . The collective operators 
are obviously rup, and the "conventional" monopole force 
Vomm = Vo{J2iTip)^- Granted that Vq would be scaled 
by A^^, as a centroid, it would have to be calculated for 
the occupied orbits only, which would have to be varied 
to detect a minumum. In other words: there is no col- 
lapse, but the operator is useless: All the information is 
contained in Vq, which is not a constant, but the result 
of HF variation from nucleus to nucleus. 

However, it is clear that W in Eq. (^^ is the normal- 
ized version of the mm operator, and that it is quite a 
useful. 



III. PAIRING MISCELLANEA 

The treatment of the monopole term is new, and I can 
only hope that the techniques and results in M will stim- 
ulate further interst in Tim- 

The quadrupole force is very old, and nobody seems to 
have much problems concerning it. Here I only note that 



the quadrupole terms involving 2Huj jumps are quite sup- 
pressed in the reahstic force: They happen to be exactly 
what is needed to produce the correct effective charges in 
perturbation theory. Therefore we are left basically with 
a q ■ q force of Elliott's type, associated with the SU(3) 
symmetry p5| ]. 

The other big multipolc terms in the collective Hamil- 
tonian are very much what they should be: octupole, 
hexadecapole and ar forces. There is also a huge 7 = 10^ 
center of mass contribution, that is nicely decoupled, as 
it should, to ensure momentum conservation. 

We are left with the pairing terms. Contrary to 
quadrupole — which seems to be naturally and satisfac- 
torily included in mean field formulations — pairing has 
to be put by hand. Since the old version is not suflicient 
for state of the art calculations, there are many proposals 
to replace it. Furthermore, there is also much interest in 
the other pairing terms. Let us start with a look at these. 

The very large F = 10 one is the ST = 10 part of a 
pairing force in LS scheme (F = 01 is indeed the ST ~ 01 
part). In jj coupling, the term is massively dominated by 
the matrix elements between the largest I ■ s partners (e. 
g., 77/2/5/2 in p = 3). In spite of its strength it is much 
suppressed by the splitting between the orbits. Further- 
more, its T = nature makes it rapidly inefficient when 
moving away of T = nuclei. It is also at the origin of a 
common misconception regarding the Wigner term, that 
deserves a digression. 

The Wigner term is the piece of the force that goes 
like T. It simply comes from the T(T + 1) symmetry 
energy which is very strong, and produces a cusp in the 
binding energies at T = 0. Therefore, before deciding 
that F = 10 (or some JO) pairing is responsible for the 
Wigner term, people are asked to check that the coef- 
ficient of the symmetry energy (a centroid of type b in 
Eq. (y)) has been kept constant. 

The F = 20 term is a puzzle. I do not know what to 
make of it. 

There are no other pairing terms that may claim a col- 
lective status, in particular, the F = 21 term does not 
seem to amount to much. I would strongly urge people 
to stop using it. 

Finally, let me return to F = 01 pairing. As men- 
tioned, there are various new candidates. Once the ones 
that do not make sense are discarded (I have in mind 
the 6 force) , my feeling is that something like the version 
that has been proposed here will be either accepted or 
independently discovered. Modifications may be neces- 
sary, because there are problems (discussed in M), but 
the proposal seems quite sound. 



Refs. [[16|-|18[ contain discussions about the influence 
of F = 01, 10 pairing on backbending rotors, that can be 
summed up as follows: 

Although the energetics of the yrast band are strongly 
affected by the pairing modifications, the other properties 
are not, since the wavef unctions change little. 

Nuclear physics is not confined to backbending ro- 
tors. In other regions, pairing will dominate the coupling 
schemes, with quadrupole acting as a pertubation. In 
general, it seems quite likely that nuclear structure will 
remain a field in which pairing and quadrupole are the 
main players, with monopole acting as referee. 



[1] At the conference I was originally supposed to talk 
about pairing, and I tried to convey the main ideas 
in a few minutes. For an expanded version of the 
main presentation, which dealt with level densities see 
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